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The eikonal, partial wave (PW) Lippmann- Schwinger, and three-dimensional 
Lippmann- Schwinger (LS3D) methods are compared for nuclear reactions that 
are relevant for space radiation applications. Numerical convergence of the 
eikonal method is readily achieved when exact formulas of the optical potential 
are used for light nuclei (A <16) and the momentum-space optical potential is 
used for heavier nuclei. The PW solution method is known to be numerically 
unstable for systems that require a large number of partial waves, and, as a 
result, the LS3D method is employed. The effect of relativistic kinematics is 
studied with the PW and LS3D methods and is compared to eikonal results. 
It is recommended that the LS3D method be used for high energy nucleon- 
nucleus reactions and nucleus-nucleus reactions at all energies because of its 
rapid numerical convergence and stability. 


1 Introduction 

The space radiation environment is composed of solar particle emissions and ions produced 
from supernovae distributed throughout the galaxy [1, 2]. Solar particle events, including 
both coronal mass ejections and solar flares, are composed of mostly protons with energies 
that can exceed several hundred MeV. Galactic cosmic rays originate from the shock waves 
of supernovae and consist of protons and heavier ions with energies that reach hundreds of 
GeV per nucleon. Radiation transport codes are used to describe the transport of ions, and 
secondary particles produced from nuclear collisions, from the space radiation environment 
through shielding materials. Space radiation transport codes require cross sections for 
the numerous nuclear reactions that occur as a result of collisions of nuclei in the space 
radiation environment with nuclei in the shield. NASA’s deterministic transport code, 
ffZETRN [3-5], currently transports all ions up to nickel — where, thereafter, incident 
particle fluxes are negligible [6] — with energies that extend from MeV to hundreds of 
GeV per nucleon through shielding materials. Efficient, accurate codes are needed for the 
computation of nuclear cross sections due to the large number of nuclear reactions that 
occur at these energies. 

The Lippmann-Schwinger (LS) equation is an expression for the scattering transition 
amplitude [7]. Scattering amplitudes can be obtained by either solving the LS equation 
or by employing some approximation, such as the eikonal method. The elastic differential 
cross section is computed from the absolute square of the scattering amplitude, and the 
total cross section is related to the imaginary part of the forward scattering amplitude. 
The elastic cross section is obtained by performing the angular integration of the elas- 
tic differential cross section, and the reaction cross section is found from the difference 
between the total and elastic cross sections. 

For elastic scattering, it is convenient to use the ground state and excited state pro- 
jectors to write an equivalent set of coupled equations for the LS equation: the elastic 
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scattering equation and the defining equation for the optical potential [8, 9]. This is use- 
ful for the study of elastic scattering because all of the complications that arise from the 
intermediate excited states are not explicitly found in the elastic scattering equation. The 
input into the elastic scattering equation is the optical potential, which can be expressed 
in an infinite series of nucleon-nucleon (NN) transition amplitudes, £nn- If the transition 
matrix is written for ground states of the projectile and target, then, in the factorization 
approximation, the first order optical potential is proportional to £nn and the nuclear 
densities of the projectile and target [10-12], The model of £nn used in the present work 
is parameterized to NN total cross sections, slope parameters, and the real to imaginary 
ratios of the transition amplitude. Nuclear charge density distributions are obtained from 
electron scattering experiments [13, 14], Matter densities of nuclei are found from nu- 
clear charge densities by factoring out the charge distribution of the proton. The internal 
charge structure of the proton is not taken into account in this analysis; instead, nucleons 
are treated as point particles. Harmonic well densities are typically used for lighter nuclei 
because of the Gaussian-like decay of the nuclear charge density as a function of radial 
distance. Wood-Saxon densities, also known as two-parameter and three-parameter Fermi 
densities, are better suited for heavier nuclei, where the nuclear charge density is relatively 
constant before decreasing to zero at larger radial distances. 

The two most common ways of solving the LS equation are to use the eikonal approx- 
imation or the method of partial wave (PW) decomposition. The eikonal approximation 
was first introduced by Moliere and systematically developed by Glauber in the treatment 
of many-body nuclear reactions with a quantum collision theory of composite objects [15]. 
The eikonal approximation can be derived by assuming high energy and small angle scat- 
tering, which leads to a linearized propagator in the LS equation from which the eikonal 
scattering wave function may be obtained [7]. The scattering amplitude is determined 
from the eikonal phase factor, which is a function of the the optical potential [5, 16-20]. 

Besides being an approximation, a drawback of the eikonal approximation is that it 
may be numerically inefficient for the evaluation of the cross sections for a given opti- 
cal potential. In the position-space representation, the optical potential, U( r), is given 
by a 6-dimensional integration for heavy ion collisions [5]. Therefore, the eikonal phase 
factor depends on a 6-dimensional integral in the position representation of the optical 
potential and an additional integration variable over a coordinate in the scattering plane. 
The numerical integration over 7-dimensions in the position space representation is inef- 
ficient when an analytic expression of the optical potential is not known. It is desirable 
to use exact formulas for the optical potential when analytic expressions of the optical 
potential can be found. The current work implements expressions of the optical potential 
for nucleon-nucleus (NA) and nucleus-nucleus (AA) scattering utilizing harmonic well nu- 
clear matter densities for light nuclei (A < 16), and the optical potential is expressed in 
momentum space for cases where no analytic expression can be found (A > 16). 

The present work expresses the optical potential as a function of momentum transfer 
and reduces the 7-dimensional integration of the eikonal phase function to 2-dimensions, 
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which significantly increases numerical efficiency. An additional advantage of this formu- 
lation of the eikonal phase function is that the momentum space representation of the 
optical potential may be used for any nuclear density and transition amplitude, provided 
the Fourier transforms may be evaluated or the explicit forms of the the densities are 
known in momentum-space. 

The LS equation may also be solved via the method of partial wave decomposition 
[7, 21], where the transition amplitude is expanded in an infinite series of functions of 
relative momenta and angular dependent spherical harmonics or Legendre polynomials. 
After projecting out the angular dependence, the transition amplitude is solved for a given 
partial wave. Once the partial wave solutions are found, the full solution for the transition 
amplitude is found by re-summing the series, which is terminated when some pre-defined 
tolerance of precision is reached. 

The PW method is known to become numerically unstable for reactions that require 
many partial waves [7], which is not only limited to high energy NA reactions (GeV/n) 
but also includes AA reactions at relatively low energy per nucleon (hundreds of MeV/n). 
The numerical instability can be traced back to highly oscillating Legendre polynomials in 
the PW expansion and large on-shell momenta for elastic reactions, where contributions 
to the transition amplitude tend to be localized. 

Although there are numerical limitations associated with the PW method, the full 
three-dimensional Lippmann-Schwinger (LS3D) solution method circumvents the neces- 
sity of using highly oscillating Legendre polynomials [22-26]. Most of the LS3D studies 
have consisted of NN interactions [22, 23, 25] with the exception of Rodriguez-Gallardo 
et al. [24] who studied NA and AA reactions at relatively low energies and Liu et al. [26] 
who studied three-body reactions. This demonstrates the validity of the method and can 
be compared to results generated with the PW method since few partial waves are needed 
for such reactions. The LS3D method is compared to the PW and eikonal methods for 
NA and AA reactions with energies extending from 150 MeV/n to 20 GeV/n. 

The eikonal solution method is a non-relativistic approximation; however, when ener- 
gies become sufficiently high, relativistic effects will be manifested in the elastic differential 
cross section. Relativistic kinematics are needed for high energy reactions and are easily 
incorporated into the momentum-space representation of the PW and LS3D equations, 
where the momentum is simply a number instead of a spatial derivative operator, as in 
the position-space representation. At relativistic energies, the PW and LS3D models will 
agree if convergence of the partial wave solution is reached, but both methods should 
differ from the eikonal results. In the low energy limit, the eikonal method should break 
down and begin to diverge from the PW and LS3D results because of the approximations 
inherent in the eikonal method. To examine the effect of kinematics, model results are 
compared for various nuclear reactions at relativistic and non-relativistic energies. 

At relativistic energies, the inner structure of the nucleons may be probed. The mul- 
tiple scattering theory (MST) upon which the model of interaction is based and the NN 
transition amplitude do not account for the inner structure of the nucleons. The compli- 
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cations associated with the inner structure of the nucleons are assumed to be included in 
the parameterizations to experimental NN transition amplitudes. 

In this paper, exact formulas are presented for the optical potential in the position 
space representation for light nuclei. A momentum-space representation of the eikonal 
approximation is presented and used for heavier nuclei. Both methods are utilized in 
the present work and are significantly more efficient than performing the usual eikonal 
calculation in the position-space representation. The PW and LS3D methods are solved 
with non-relativistic and relativistic kinematics, and comparisons of the models are made 
for reactions that are relevant to space radiation. Additional comparisons to experimental 
data are presented. Based on the results presented herein, it is recommended that the 
LS3D method be used for high energy NA reactions and AA reactions at all energies 
because of its rapid numerical convergence and stability. The effect of the kinematics for 
projectiles and targets with equal masses and extensive comparisons to experimental data 
will be communicated in subsequent manuscripts. 

The present work is organized as follows. In section 2, a theoretical overview of the 
LS equation, MST, the elastic scattering equation, and the optical potential are reviewed. 
This is followed by a discussion of the eikonal, PW, and LS3D methods. The numerical 
solution methods are discussed in section 3, and model comparisons are given in section 
4. The conclusions are stated in section 5. 


2 Theory 

2.1 Lippmann-Schwinger Equation 

The LS equation is an expression for the scattering transition operator — the fundamental 
quantity that is used to evaluate the elastic differential, elastic, reaction, and total cross 
sections for nuclear reactions — and is given as 

T —V + VG+T, (1) 

where V is the sum of residual two-body interactions (potentials) for the projectile-target 
system, and Gq is the unperturbed two-body propagator [27]. In momentum-space, the 
transition matrix element may be written 

<k'|T|k) = (k'|V|k) + (k'VGJTIk), (2) 

where k and k' are the initial and final relative momenta of the projectile-target system 
in the center of momentum (CM) frame, respectively. Using the notation T(k',k) = 
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(k'|Xjk), the momentum-space integral equation representation of the LS equation is 

T( k', k) = V(k', k) + J V(k', k ")G+{k", k)T( k", k)dk", (3) 

where /c = |k| and completeness over momenta has been used, 

i = J dk"|k")(k"|. (4) 

For elastic scattering, the scattering amplitude, /($), is related to the transition matrix 

by [7] 

m = — ( 5 ) 

where k — |k|, 6 is the scattering angle, the density of states, p, is given by 

p = k 2 dk/dE (6) 

and E is the energy. For non- relativistic (NR) kinematics, E = k 2 / 2 p, where p = 
( mprriT)/(mp + mp) is the reduced mass, mp is the mass of the projectile, and m T is 
the mass of the target. When using relativistic (REL) kinematics, E = \Jk 2 + m 2 p + 
y/k 2 + rrij,. 

The total cross section is proportional to the imaginary part of the forward scattering 
amplitude, cr to t = (47r/fc)Im [f(9 = 0)], and the elastic differential cross section is found 
from taking the absolute square of the scattering amplitude da/dfl = \f(6)\ 2 . One simply 
integrates the differential cross section over the solid angle to obtain the elastic cross 
section, cr el = f(dcr/dfl)d{l. The reaction cross section, <r re , can be found by subtracting 
the elastic cross section from the total cross section, a re = cr to t — ow 

2.2 Multiple Scattering Theory 

The fundamental theory which describes the interaction is MST, which was developed 
because the Born series proved to be inadequate due to the strong nature of the nuclear 
force. Watson showed that the unperturbed Hamiltonian can be separated from the sum 
of nucleon-nucleon interactions (residual interaction) between the projectile and target 
such that the leading term in the Born series can be expressed as the sum of pseudo two- 
body Watson-r operators [28] . For AA scattering, the transition amplitude from equation 
(1) is written [27, 29] 

Ap Ap 

t = p) 

*= 1 j = 1 
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where 


Ap Ap 

Tij = + Tij G+ 'y ^ 'y^ Tki ■ 

kjti Ijtj 


( 8 ) 


Ap and Ap are the number of nucleons in the projectile and target, and the Watson-r 
operator is 

Tij = % + v ij G^r ij . (9) 

The objective of the MST is to express the transition amplitude in terms of the free 
NN transition amplitude, = i\j + v^gt ^ , where g is the free two body propagator. To 
that end, the pseudo two-body operator is written [27, 29] 


Ti 




— tij + tij(G 0 g)Tij. 


( 10 ) 


The impulse approximation is often used and is based on the idea that the projectile 
has sufficiently high energy such that the target nucleon is treated as if it were free; that 
is, Tij ps tij , which leads to the following transition amplitude [27, 29]: 


A p Ap A p Ap A p Ap 

-^impulse tij + tijG 0 tkl + m ■ ■ 

i = 1 j = 1 i = 1 j= 1 k^i l^j 


( 11 ) 


The propagator Gq includes the ground states and all intermediate excited states of the 
projectile and target. As a result, equation (11) is difficult to solve. 

If elastic differential cross sections are sought, then it is convenient to separate the 
LS equation into an equivalent set of coupled equations with the ground state (P) and 
excited state (Q) projectors [8, 9]. The elastic scattering equation is 


T = U + UPG^PT, 
and the optical potential is given by 

U = V + VQG+QU. 


( 12 ) 


(13) 


All of the intermediate excited states have been included in the optical potential. With 
that in mind, the Watson series for the optical potential may be written in a form similar 
to equations (7) and (8) [29], 

A p Ap 

U = '£'£u ij , (14) 

*= 1 3 = 1 

and 

A p Ap 

Uij = fij + f l3 QGy q y y u kh (is) 

k^i l^j 
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with 


(16) 


T{j — V%j QGq Ql~ij • 

Again, f may be expressed in terms of the free two-body amplitude [27, 29] 

7~ij tij tij (QG q Q g)T %J , 


(17) 


where the tilde notation, f, indicates that the pseudo two-body operators have been 
modified by the excited state projector, Q. If the impulse approximation is used (f^ ~ t l3 ) , 
the Watson series for the optical potential becomes [29] 

Ap Ap Ap Ap Ap Ap 

^impulse k,QG+Q EE tki + • • • • (18) 

i= 1 3 = 1 i=l J=1 


Finally, if the higher order terms are neglected, and only single scattering is considered 
[27, 29]: 

Ap Ap 
i = 1 t=l 


r/ lst 

impulse 


At higher energies, it is assumed that higher order terms do not contribute significantly 
more than the single scattering term. 

In all of the results that follow, U = ^jmp ulse , an d the elastic scattering equation 
becomes 


T(k',k) = t/(k',k) 


P'(k / ,k // )T(k /, ,k) 
E{k) — E(k") + ie 


( 20 ) 


where E is the energy and ie ensures outward scattering boundary conditions. 


2.3 Matrix Element of Optical Potential 

The explicit form of the optical potential must be defined before the solution of the LS 
equation can be found. The matrix element for the optical potential is written [10] 

A p Ap 

t/(k',k) = (k'-,^^\U^^;k) = ^^(k'-X P K T \t,AK P ^ T -,k), (21) 

i = 1 3 = 1 

where \4>q P ) and \4>q T ) are the ground states of the projectile and target, respectively. 
Summing over individual NN interactions, the matrix element becomes 

e(k'.k) = A P A T (k'-,^-^\t\4>i P ^ T -,k). (22) 

Single scattering between a projectile and target nucleus is shown in Fig. 1. k and k' 
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are, respectively, the initial and final momenta of the projectile nucleus in the AA CM 
frame. Likewise, — k and — k' are the initial and final momenta of the target nucleus in 
the AA CM frame. The internal virtual nucleon momenta over which the integration is 
performed are denoted pi and P 2 . P p and P p are the momenta of the core of the Ap — 1 
virtual projectile and Ap — 1 virtual target nuclei, respectively. Inserting completeness 
over the initial and final internal nucleon momenta and using conservation of momenta 
in equation (22) gives 

U (k', k) = ApAp J (k , ;0o T 0o P |p , 2 pj)(p2pj|t|p 2 pi)(p 2 pi|k;0o T 0o P ) ( 23 ) 

x dpjdp' 2 dpic/p 2 <5(pj - Pi - q)5(p 2 - p 2 + q), 

where q = k' — k is the momentum transfer. 

After evaluating the delta functions, equation (23) becomes, 

f/(k',k) = A P A T /(^|p 2 -q)(p 2 |^)(^|pi + q)( Pl |^) (24) 

X (k'; p 2 - q, Pi + q|t|k; p 2 , pi)fip,dp 2 . 

The notation above includes initial and final relative momenta of the projectile/target 
system in the CM frame, k and k', and the internal nucleon momenta relative to the CM 
of the (projectile or target) nucleus, denoted p,, with i = 1,2. These momenta (p«) are 
associated with the relative kinetic energy of the nucleons within the projectile (or target). 
The semicolon is used to show that the relative internal momenta are in a different space 
than the relative momenta. In the discussion that follows, k and k' do not change and 
are included for the convenience of later manipulation. 

2.3.1 Optimum Factorization 

Next, the optimum factorization for NA reactions in the CM frame [12, 30-32] is extended 
to AA reactions in the CM frame and is shown in Fig. 2, where the symmetric coordinates 
have been chosen for time reversal invariance [12, 30] and to allow for easier calculations. 
The initial momentum of the projectile nucleon is pi — q/2, and the final momentum is 
Px + q/2. Similarly, the initial momentum of the struck target nucleon is p 2 + q/2, and 
the final momentum is p 2 — q/2. It is important to note that the momenta were chosen 
such that q = k' — k and is the reason why the signs are reversed for the momentum 
transfer of the target nucleon. 

The wave functions will eventually be related to the nuclear charge densities. Thus, 
the arguments of the wave functions will be written as arguments of the nucleon momenta 
relative to the cores of the Ap — 1 and Ap — 1 nuclei. Let P 1 be the momentum of the 
projectile nucleon relative to the core of the A P — 1 nucleus; likewise, P 1 is used for the 
target nucleons. The initial momentum is denoted i and the final momentum is /. In 


addition, if K = (k + k')/2 and momentum conservation from Fig. 2 are used, then the 
core momentum of the A p — 1 virtual projectile nucleus can be expressed as Pp = K — pi. 
It also follows that Pp = — K — P 2 . 

The relative momentum between particles a and b with momenta p a and p & and masses 
m a and mj, is defined "P = (m a pb — rribp a ) / (m a + mb). The initial relative momentum of 
the projectile nucleon can be found by inspection of Fig. 2. Let p a be the momentum 
of the core of the Ap — 1 nucleus (p a = Pp), and let p& be the initial momentum of the 
projectile nucleon (pf, = pi — q/2) with m a = Ap — 1 and = 1. This leads to an initial 
relative momentum of the projectile nucleon relative to core of the Ap — 1 nucleus of 



Similarly, the other relative momenta may be found, 



= Pi + 



K 

A? 



and 



(25) 


(26) 


(27) 


( 28 ) 
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The optical potential becomes 


t/(k',k) = A p A t J MV T t ,Vj)pp(T^-,Vl 


(29) 


/i / q q, ,, q q. 

X (k ;p 2 - -,Pi + ~\t |k;p 2 + ^Pi “ -)dpidp 2 , 


where p is the one-nucleon (proton or neutron) density matrix of the nucleus and is related 
to the nuclear charge density, p, by [12, 30] 


p(q) = P I 


A - 1 


A 


~Q = P y - 


A-l 

A 


q . 

2’ y + 


A-l 

A 


b dy - 


(30) 


where y is a dummy integration variable. 

It is convenient to make a change of variables by defining the following: Pi = pi — 
K f Ap and P 2 = p 2 + K /At- Then the matrix element of the optical potential becomes 


U (k', k) — A P A T J p T |^P 2 — 
x Pp f Pi + 


Ai 


An 


q p , 

2 ,P2 + 


Ai 


An 


5 |rfp,dp 2 


(31) 


Ap — 1 

q P, 

1 

l 

Ap 

2 ,P 

Ap 


q 

. 2 

„ , „ q K ^ q K , „ „ q K q K. 

x <k ;P2 ~ 2 “ t t ' Vi + 2 + i ; |t|k;P2 + 2 " i? Pi ” 2 + 

where the transition amplitude is evaluated in the AA CM frame. The transition ampli- 
tude in equation (31) can be written 


„ q K ^ q K . „ q K ^ q K. 

(k ; P 2 — — Pi + — + — |f|k; P 2 + — — , P 1 — — + — ) 


Ai 


A P r '“’ ‘ ' 2 At a 2 
= (k^k^k'Jtlk, k 2 ;ki>, 


Ai 


(32) 
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where the initial of final projectile and target nucleon momenta in the AA CM frame are 


q K 

P ‘~2 + T p 

(33) 

„ q K 

P2 + 2 " T r 

(34) 

ei K 

Pl + 2 + T P 

(35) 

P q K 

2 2 A t 

(36) 


The current work uses parameterizations of NN amplitudes to predict NA and AA 
differential cross sections. Consequently, the quantity that allows for the transformation 
of the transition amplitude from the AA CM frame the NN CM frame, known as the 
Moller factor ( 77 ) [7, 33], must be found. 

The Moller transformation factor can be found by examining the Lorentz invariant 
(LI) cross section, which is written [7, 12] 

da = J '^r'^rKE ~ E')5 ( P - P')^|M 12 | 2 , (37) 

with 

M 12 = \J E[E' 2 W EiE 2 , (38) 

where E and P are the initial total energy and momentum of the system, E' and P' are 
the final total energy and momentum of the system, Ei and E 2 are the total energies of 
the two particles before the collision, E[ and E' 2 are the total energies of the two particles 
after the collision, ki and k 2 are the initial momenta of the two particles, kj and k ' 2 are 
the final momenta of the two particles, and B = \ a / A(s, mf, m|), where 

A(s, m\ , ml) = (s — m\ — ml) 2 — Am\ml. (39) 

Since, dlti/Ei, 5(E — E')5(P — P'), B, and da are LI, then M 12 is LI. Thus, 

[Mi 2 \ aa = [M i2 ] NN . (40) 

The initial and final relative momenta of projectile and target nucleons in the NN CM 
frame can be expressed as 


k = 


1 

2 




(41) 
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and 


k'= -- 


P,-5-K 


Pi+ 2 + Jp 


2 At 

Using equation (38), the AA CM frame is related to the NN CM frame by 
v /£ N (k' 2 )£ N (k' 1 )<k , ; k' , ki|t|k; k 2 , k 1 ) AAV /U N (k 2 )U N (k 1 ) 


= y -S N (k / )-S N (— k')(k / ; k', -k'|t|k; k, -k ) NN ^ E N (k)E N (~k) 
Therefore, < t >aa= V < t >nn, where 


V = 


1 E N (k')E N (-k')E N (k)E N (-k) 


U N (k 2 / )A N (k 1 ') J UN(k 2 )U N (k 1 )’ 

The notation can be simplified by defining 


Y = -(k + k') = - 
2 V ; 2 


K ^ K 

^ -P2 + a; + Pi 


such that 


and 


(k'; k 2 , k[\t\k] k 2 , k x ) AA = ?7t NN (q; -2Y, 2Y). 


U(k',k) = A p A t I p T (P 2 

Ap — 1 


Ap — 1 


x Pp\ P i + 


A, 


Ai 


S p 

2 ,Pl 


q p , 

j,P2 + 


At — 1 


Ai 


q 

2 


Ap — 1 

Ap 


(42) 


(43) 


(44) 


(45) 


(46) 


(47) 


2 }^NN(q; -2Y,2Y)dPic/P 2 . 


The NN transition amplitude may be written as a function of q and YfPnPo) without 
loss of generality. 

Since the size of the nucleus is much larger than the range of the NN interaction, the 
transition amplitude varies more slowly with Px and P 2 when compared to the one nucleon 
density matrix for the nucleus [12, 30]. For this reason, the NN transition amplitude is 
evaluated at fixed values of Pi and P 2 , namely Pi = Pio and P 2 = P 2 o- The optimum 
values of Pio and P 2 o are estimated by expanding the transition amplitude in a Taylor 
series and requiring second-order terms (and higher) to be zero. After the optimum values 
are found, the optical potential for elastic scattering is stated for central potentials. 
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The transition amplitude may be expanded as 


^NN(q; 2Y) « t N x(q; — P20 + - 7 — b Pio) ( 48 ) 

Zip Zi p 

+ 0T: , ^NN(q; -r — P20 + -r — b P10) (P2 — P20) 

ar^20 /it /ip 

+ — %N(q; -7 P20 + ~7 b Pio)(Pi — P10) + ^(Pi> P2), 

O Jr iQ /±p /ip 

where 6 represents higher order terms. Using equation (30) and the Taylor series expansion 
leads to the following matrix element of the optical potential: 


f/(k',k) = A P A T rjt NN (q;, - P 2 o + yr 1 + Pio )PT(q)pp(q) 

Zip Zip 

' A d . K K „ " 

A P A T „ ^ NN (q; -7 P20 + —7 b P 10) 

Or 20 /It /ip 


(49) 


X 


P 20 )Pt(P 2 , q)pp(Pi, q)dPidP 2 


+ 


' , d . K K ^ 

A P A T „ t NN (q; — P20 + -7 b Pio) 

Or 10 Zip Zip 


X 


(Pi - Pio)pt(P 2 , q)p P (Pi, q)dPidP 2 


+ . . . . 


The optimum values for P 2 o and Pio can be estimated by requiring the second order 
terms in the Taylor series expansion to be zero [12]; that is, 


(P 2 - P2o)/MP2,q)/5p(Pi,q)<iPi<IP2 - 0 


( 50 ) 


and 


/(Pi - Pio)/ 5 r(P 2 ,q)pp(Pi,q)dPidP 2 = 0 . 

Equation (50) can be written as 


(51) 


P2oPr(<?)Pp(<?) - / Pp(g)pr(P 2 , q)P 2 'lP2- (52) 

The single particle density matrix for the nucleus is derived from the ground state wave 
function, which is a symmetric function of its coordinates and invariant under the parity 
transformation (P 2 — >■ — P 2 ) [12], since the strong interaction preserves parity. Likewise, 
the entire integral is invariant under a parity transformation and evaluates to zero; thus, 
P 20 — 0. The same argument follows for the second term of the matrix element of the 
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optical potential and Pi 0 = 0. Using these values for P 10 and P 2 o leads to 


U (k', k) — A P A T r]t NN 


q; 


A P + Ap 

ApAp 


K 


PtMppM- 


( 53 ) 


Only central potentials have been considered; that is, the effects of spin are ignored. 
As a result, the interaction must be a function of the possible scalar products in equation 
(53): q - q, q - K, and K ■ K. Using |k| = |k'| for elastic scattering, one can show that 

£/(k',k) = ApAp7]tis!Ki(eNN, q)pt(q)pp(q), (54) 


where cnn is the NN CM kinetic energy. The current work uses an NN transition am- 
plitude that is evaluated at the energy available per nucleon that is related to the fixed 
beam energy [11]. 

In the next section, the eikonal approximation is discussed, and exact formulas for the 
optical potential are given in position-space with the transition amplitude and harmonic 
well densities described in Appendices A and B. 


2.4 Eikonal Solution Method 

The eikonal approximation is used for high energy, small angle scattering to calculate 
elastic, reaction, total, and elastic differential cross sections [7, 34], To compute cross 
sections with the eikonal method, one solves for the eikonal scattering amplitude, /(#), 
which is given as [7] 


f(6) = j J J 0 (2ksm(9/2)) [e ix ^ b) - l] b db : (55) 

o 

where k is the relative momentum of the projectile-target system in the CM frame, Jo is 
the ordinary cylindrical Bessel function, 6 is the scattering angle in the CM frame, b is 
the impact parameter, and y(/c, b) is the eikonal phase shift function, which is obtained 
by integrating over the optical potential, U(b, z ) [7]: 

OO 

x(k,b) = - ^ J U{b,z)dz. (56) 

— OO 

The ^-integration is taken to be in the same direction as the initial wave vector of the 
incident projectile. The elastic differential cross section is found by taking the absolute 
square of the scattering amplitude, da / dVt = \f(9)\ 2 . 

Since the eikonal approximation satisfies the optical theorem, the total cross section 
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is proportional to the imaginary part of the forward scattering amplitude, 


47T [ 

o tot = f(0 = 0) = 4vr J [ 1 

o 


- e- Im * cos (Re x)}b db. 


( 57 ) 


The total elastic cross section, cr e i, is obtained by integrating the elastic differential cross 
section [7], 


[1 — e Imx cos(Rey)]5 db 
o 

OO 

- 2vr /[l- e _2Imx ]6 db , 
o 

and the reaction cross section, <r re , is found by using cr re = a to t — o e p thus, 

OO 

<r re = 27 t J[ 1 — e _2Imx ] 6 db. (58) 

o 

All of the cross sections described above are functions of the eikonal phase function 
and depend on the optical potential, as shown in equation (56). For AA scattering, the 
optical potential may be expressed as [5, 16] 

U(r) = A P A T J t NN (|r NN |K(M)MM) ^dr NN , (59) 

where A is the number of nucleons, P represents the projectile, T represents the target, 
£nn is the NN transition amplitude, and p is the nuclear matter density, tnn is the vector 
between a nucleon in the projectile and a nucleon in the target, vp is the vector that 
extends from the center of the projectile nucleus to a nucleon in the projectile, iy is 
the vector between the center of the target nucleus to a nucleon in the target, r is the 
relative distance between the centers of the projectile and target nuclei, R = r + rp is 
the distance from the center of the projectile to a nucleon in the target. The distance 
from the center of the projectile nucleus to a nucleon in the projectile may be expressed 
as r P = r + R = r + r T + r NN , which, when substituted into equation (59), leads to 

U(r)=ApA T J tNN(| r NN|)Pp(| r + r T + r NN|)PT(| r T|)^ r T^ r NN- (60) 

(Note that |r| = \fb 2 + z 2 in the cylindrical coordinate system.) The vectors used in this 


= l^m dSl = 41 
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discussion of the AA optical potential are illustrated in Fig. 3. 

In the next two sections, analytic expressions for the optical potential are presented 
for nuclei with A < 16, where harmonic-well nuclear matter densities have been used. 
Harmonic well densities [5] are given by 



(a + f3r 2 ) exp 



and the NN transition amplitude is [5] 


tNN(r) = rex p 


25(e) 


(61) 


(62) 


(See Appendices A and B for a full description of the densities and transition amplitudes 
used in the current work.) 


2.4.1 Nucleon-Nucleus Optical Potential 


For NA collisions, the single projectile nucleon charge density is taken as a Dirac delta 
function, and the harmonic- well nuclear matter density is used for the target. Thus, the 
exact optical potential from equation (60) may be expressed in the following form with 
A P = 1: 

U (r) = (C 0 + Cir 2 ) exp [~C 2 r 2 ] (63) 


with 


and 


where 


3/2 


Cn — tAi 


7 T 


H 1 + IT 2 


OCT + 


3 Pt 

2(/4l + /ha). 


Ci = 


r/3 r A r /i 2 7 r 3 / 2 


(/H + /D ) 7 / 2 ’ 
C2 = pi — 


f4 


(pi + p 2 ) ’ 


pi = ^ and /i 2 = 


4s|' 


(64) 

(65) 

( 66 ) 

(67) 


2.4.2 Nucleus-Nucleus Optical Potential 

The calculation of the optical potential is repeated for AA collisions with harmonic- 
well nuclear charge densities for both the projectile and the target, which results in the 
following formula for the optical potential: 
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where 


U{r) = (A 0 + Air 2 + A 2 r 4 ) exp[-A 3 r 2 ], 


( 68 ) 


Ao 


27t N ApAj' 

e 



X 


otT Ai 3 . N 

— 2 ^ 40 + A 2 oip ) 

15A 2 /3t 
" 80 2 


A 


2ttN ApAp 
0 



x 


°^ + ^(3A 2 /3 T - Aa T A 2 S) 


+ 7 ^ (PtAi5 + A 2 ap8 — 5f3pA 2 ) 

58 2 A 2 /3t 
+ 20 3 ’ 


and 


where 


2ttN A pAp f 7r 

A2 = ? V» 

~A 2 /3 T 8 2 (3 t A 2 8 3 A 2 /3 t 8 4 ' 

X 20 2 + 20 4 ’ 



AT 


2ttt 

K 



K = 


1 

44 


1 

25’ 


e = 


1 

Asp 


+ 8, 


(69) 


(70) 


(71) 


(72) 


(73) 

(74) 

(75) 
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5 


1 


1 


(76) 


and 


4 s 2 p 16 SpK 1 


Ai 



3 (3p 

~1k' 


(77) 


Pp Pp Pp 

2 32K 2 Sp 4 KSp 


(78) 


These results have been derived independently by Townsend and Bidasaria [35, 36]. In 
the next section, the optical potential is written as a function of momentum transfer, 
which leads to more efficient evaluation of cross sections when exact expressions for the 
optical potential are not known. 


2.4.3 Optical Potential in Momentum-Space 

In section 2.2, it was shown that the phase function is obtained by integrating the optical 
potential in the position-space representation. In this section, the phase shift function is 
written such that the optical potential is expressed as a function of momentum transfer 
q. By doing so, the number of integration dimensions will be significantly reduced. 

To begin, the NN transition amplitude and nuclear charge densities in equation (60) 
are replaced with their Fourier transforms 


U(r) 


A pAj 1 

(27 r) 6 


dr T dr NN dqi dq 2 dq 3 [t NN ( | qi | ) Pt ( | q 2 1 ) Pp ( | q3 1 ) 


x e 


-*qi r N N e -*q 2 r T 


e -jq3-(r+r T +r N N)j 


(79) 


where the Fourier transforms are given by 

^ NN ( r ) = Ay J t N N(q)e _ * q ' r dq, (80) 

and 

P( r ) = J p(q)e" 4 q ' r dq. 

The normalizations for the Fourier transforms of p{ r) and 
traditional normalization of nuclear charge densities. Next, 
is performed with the Dirac delta function, 

■5 (a) =(2 < 82 > 


(81) 

^NN( r ) differ because of the 
integration over and tnn 
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which results in 


U{y) = A p A t / t N N(|qi|)M|q2|)pp(|q3|)6(qi - q 3 )5(qi + q 2 )e * q3 ' r dq x dq 2 dq 3 . (83) 


/ 


After evaluating the delta functions, the optical potential is reduced to integration over 
the momentum transfer, 


where f/(q) = A P A r t NN (|q|)p P (|q|)p T (|q|). 

Next, the Fourier transform of the optical potential (84) is substituted into the ex- 
pression for the eikonal phase shift function from equation (56), 


where r = b+z in cylindrical coordinates has been used. The integration of 2 is performed, 
which results in 


The advantage of equation (88) is that the optical potential is in the momentum-space 
representation, and the ^-integration need not be performed. Instead, the 7-dimensional 
integral for y has been reduced to 2-dimensions over the magnitude of the momentum 
transfer, q, and the angle, 0, between the momentum transfer and the impact parameter. 
This result significantly increases the efficiency for the numerical evaluation of y. 



(84) 



(85) 


— oo 


oo 



— oo 



( 86 ) 


where the following delta function has been used, 


OO 



(87) 


— OO 


Using dq = q 2 dq sin 9 d9d(j) and evaluating the delta functions leads to the final form for 
the eikonal phase function, 



( 88 ) 


0 


0 
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2.5 Lippmann-Schwinger Partial Wave Solution Method 

The LS equation is often solved with the method of partial waves, a well-known method 
that is described in standard texts [7, 37]. An overview of the method is presented here. 

Equation (20) is solved with the partial wave method using the optical potential from 
equation (54). Any function, Z, can be expanded in a complete orthonormal set of basis 
functions, 


OO l 

z( q) = z( k', k) = y . E Z ^ k '" m m '(k')Yr(k) m 

1=0 m=—l 

°° 2/ + 1 

= '£- z -Z l (k!,k)Pi(x), (90) 

1=0 


where q = k' — k is the momentum transfer, Y) m *(P) and Y] m (fc) are spherical harmonics, 
k and k' represent the polar and azimuthal angles associated with k and k', respectively; 
Pi(x) are the Legendre polynomials, and x = cos (9 kk >), where 9 kk t is the angle between 
k and k'. After expanding T and U, one integrates out the angular dependence, which 
leads to the PW LS equation 


T t {k\k) 


U l (k , ,k) + 


OO 



0 


U l {k!,k' , )Ti{kf , ,k) 
E(k) - E{k")+ie 


k" 2 dk", 


(91) 


where 


i 


Ui(k',k) = 2vr 


U(q)Pi(x)dx. 


-l 


(92) 


A common way of solving the above equation is through the A-matrix, which consists 
of a principal value integral. When expanded in partial waves, Ri has a form similar to 
Ti above, 


Ri{k\k) 


U l (k',h) + P 


OO 



0 


Ihik^k'ORiik'^k) 

E(k ) - E(k") 


k" 2 dk ", 


(93) 


where P represents the principal value integral. The ffeitler equation may be used to 
obtain T)(/c', k), 


Tiik'^k) 


Ri{k',k) 

1 + inpRi(k'i k ) ’ 


(94) 
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where p is given in equation (6). The final step is to re-sum the infinite series, 

°° 91 - l-l 

T( q) = (k'|T|k) = E k)P,(x). (95) 

1=0 

It was recently demonstrated that elastic scattering finite sum expressions converge 
more quickly than the infinite series above [21]. The derivation is based on the idea that 
the Born approximation becomes sufficient for some large angular momentum; that is, for 
l > / max , T[(k, k ) ~ Ui(k, k). Using this idea, the finite sum expression for the scattering 
amplitude is [10, 12, 21, 38] 


^max 07 I 1 

m = E UU k) - u,(k, k)]P,(x ) - (2tt) 2 R( 9 ). (96) 

K 1=0 

The total cross section is proportional to the imaginary part of the forward scattering 
amplitude [21], 


^max 



k) - U,(k, fc)] + o*T, 


(97) 


where u^° rn is the total cross section in the Born approximation. The elastic cross section 
is derived from taking the absolute square of the elastic differential cross section [21], 




(2tt)V 
k 2 


^max 

E 


. 1=0 


2Z + 1 
47T 




+ CT 


Born 
el 5 


(98) 


where 

C"= I \ - (2*)^U(q)\ 2 i(l, (99) 

and the reaction cross section is obtained by subtracting the clastic cross section from the 
total cross section. 


2.6 Lippmann-Schwinger 3D Solution Method 

The LS3D solution method avoids the numerical difficulties associated with the PW 
method and has been used for relatively low energy reactions [22-26]. This section re- 
views the LS3D method as it has been described in the work of Elster et al. [22] and 
Rodriquez-Gallardo et al. [24], 

If one considers only central potentials in equation (3), then both T and V are scalar 
functions; that is, /(k',k) = f(k',k,k' ■ k ) for some function /, where k {k') represents 
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the angular unit vectors associated with k (k 7 ). The possible scalar products of the LS 
equation are as follows [22, 24]: 

x' = k'-k (100) 

x" = k" ■ k 
y = k" ■ k'. 

The incoming momentum, k, is taken to be in the direction of the £-axis, and the azimuthal 
angle between k and k 7 is set to zero: 0 7 = 0; therefore, y may be expressed as a function 
of x' : x", and 0" [22, 24], 


y 


= x x" + \/l — x’ 2 \/l — x" 2 COS 4 >" . 


( 101 ) 


Note that only U(k , ,k",y ) is dependent on 0"; therefore [22], 


2tt 


U{k’,x\k\x")= / U(k',k",y)d4 /', 


( 102 ) 


and the LS equation becomes [22] 

T(k’,k,x') = U(k',k,x') 


OO 1 _ 

+ fk'« dk " [ 


o -i 

The propagator may be expressed as [37] 
lim ^ ~ 


E(k ) — E{k") + ie 


v^o E(k) - E(k") +ie \ E(k) - E(k ") 


(103) 


m 6(E(k) - E(k")), (104) 


where P is the principal value. The LS equation becomes (after suppressing the on-shell 
momentum in T and V) 


T(k f , x') = U(k',x') + P 


U(k', x', k", x")T(k", x")k" 2 dk"dx" 


E{k) - E(k") 

in [ U(k',x',k",x")T(k”,x")k” 2 dk"dx"5(E(k) - E{k")). 


(105) 


Next, the delta function is evaluated, and the kinematic factors are specified. The NR 
energy is E = fc 2 /2/i, where y is the reduced mass of the projectile-target system, and 
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k is the relative momentum of the projectile-target system in the AA CM frame. The 
REL energy is E = Ep + At, where Ep = \Jm p + k 2 and At = \/rn\ + k 2 . P denotes 
projectile, and T denotes target. The delta function may be expressed as 

S(E(k) - E(k")) = ^ ~ k " ] = - k"), (106) 

\dk\(k"=k) K 


where negative momenta have been disregarded, kin represents non-relativistic (NR) or 
relativistic (REL) kinematics, and 


/i 


kin 


mpra-T 

mp+mT 


for non-relativistic kinematics 


Ep Ep 
Ep-\-Ep 


for relativistic kinematics. 


(107) 


After evaluating the delta function, equation (105) becomes 


T{k',x’) = U(k',x') + P 


u{k',x'x,x")nk " )X ") 2 

E(k) - E(k") 


(108) 


in/i 111 / U{k ' , x ' , k, x")T(k, x")kdx" . 


Equation (108) is solved with the numerical methods described in the next section. 


3 Numerical Methods 

Thus far, three different solution methods for the LS equation have been presented, in- 
cluding the eikonal appproxiation, the PW, and LS3D solution methods. This section 
describes the numerical techniques used to implement the three methods. 

3.1 Eikonal 

Exact formulas for the optical potentials in equations (63) and (68) are used for light nu- 
clei reactions (A < 16). For reactions with heavier nuclei (A > 16), the momentum-space 
representation of the eikonal phase function described in equation (88) is implemented. Al- 
though the momentum space method for the eikonal phase function is much more efficient 
than the position space calculation, additional interpolation over the impact parameter 
and momentum transfer was performed for additional numerical efficiency. Convergence 
of the total elastic cross sections was used to establish the number of Gaussian quadrature 
points used for integration. The number of Gaussian quadrature points for the eikonal 
solution method was increased to a maximum of 100 points such that the total elastic 
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cross sections changed less than 1% for all reactions with energies up to 100 GeV/n. 


3.2 Partial Wave Lippmann- Schwinger 


As discussed in the last section, one may first solve for the 7th PW of the /?- matrix, 
equation (93), and then find 7] from the Heitler equation. The transition matrix is found 
by performing the sum in equation (95). To solve for Ri, an even number of Gaussian 
quadrature points, N, is used to perform the integration of the momentum; thus, equation 
(93) becomes 

N 

R l (k',k) = U l (k',k)+J2 

k "= i 


U l (k',k")R l (k",k)k"w k „ 
E(k) - E(k") 


(109) 


where k" now represents a discrete Gaussian quadrature transformed index, and w k " is 
the corresponding weight for the integrated variable. Note that a principal value singu- 
larity will occur in the propagator when k” = k and may be removed by implementing 
a subtraction scheme [37] or the method of Sloan [39]. Sloan’s method requires that the 
pole be placed at the center of an even number of Gaussian quadrature points when in- 
tegrated. Since the pole occurs at k" = k, the limits of integration are divided into two 
parts: [0, oo) — >■ [0,2 k\ and [2k, oo). By placing the pole in the center of the interval 
[0,2k], and by using an even number of Gaussian quadrature points, the principal value 
singularity is removed. 

To isolate Ri, the above equation is written 


where 


N 


A(k',k")Ri(k' 

k "= i 

\k) = U l (k',k), 

(110) 

A{k', k") — 8 k ', k " — 

Ul(k', k")k" 2 W k n 

(111) 

E{k)-E{k") ' 


If the k' index is allowed to span the same discrete indices as k" , then the above matrix 
equation can be solved, and Ri(k",x") can be found. Observe, however, that k" are 
transformed integration points that are not necessarily equal to the on-shell momentum, 
k. Since the main focus of the current work is elastic scattering, Ri{k,k ) is ultimately 
sought. To accomplish this, an additional point is added to the sum for the on-shell 
momentum, which is not an integration point: k" = k with weight w k " = 0 for the N + 1 
point. Consequently, the matrix equation becomes 


N + 1 

Y J Kk\k")R l (k",k) = U l (k',k), 

k "= i 


( 112 ) 
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and the solution of Ri(k, k ) can be found. 

The number of partial waves needed for an acceptable tolerance of convergence is not 
known a priori. Partial waves were generated until |T) — Ui\ < 10 _4 %. 


3.3 3D Lippmann-Schwinger 

The numerical implementation of the LS3D method proceeds in the same manner as the 
PW method, but there are now two additional integration variables over azimuthal and 
polar angles. The azimuthal dependence only occurs in the potential and is integrated 
with 80 Gaussian quadrature points. As was seen with the PW method, the principal 
value integral over momenta is handled with Sloan’s method [39], and an additional on- 
shell momentum is introduced such that k" = k with wy = 0 for the N + 1 Gaussian 
point; thus, equation (108) becomes, 


T(k',x') = U(k',x')+ E 


N + 1 N TTlW U" ~/A JL//2 


k "= 1 x"=l 
N+l N 


U(k', x', k", x")T(k", x") k"*w k nw x n 
E(k) - E(k") 


(113) 


Z7r/x km U(k', x' , k" , x")T(k", x")k"w x "5k,k"- 

k"=l x"=l 


Equation (113) can be written as a matrix equation, 


N+l N 

J2 E x ") T ( k "i x ") = u ( k> i ®')» ( n4 ) 

k "= 1 x"=\ 


with 


A (k',x' : k” : x") = 


U (k 1 , x 1 , k", x")k" 2 Wk"W x n 


$k" ,k'$x" ,x' 


E(k ) - E(k") 

+ i7Tjj km U (. k x r , k" , x")k"w x "8k,k” 


(115) 


In order to solve the above matrix equation, k' is allowed to span the same discrete values 
of k"\ likewise, x' uses the same discrete points as x" . 

The system of equations can be solved to obtain T(k", x"), where the solution is given 
for the on-shell and off-shell momenta at the Gaussian transformed angles, which are 
not evenly spaced. The transition matrix can be evaluated at any angle by substituting 
T(k",x ") into equation (113). 

In order to solve equation (113), the summing indices are reduced to single indices, 
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(k',x') — > i and ( k",x ") j, by using the following relations 

i = x' + N{k' - 1) (116) 

and 

j =x" + N(k" -1), (117) 

where x', x", k f , k" represent discrete Gaussian transformed integration points (not con- 
tinuous variables). N integration points have been used for the angular integration and 
N + 1 points have been used for the integration over momenta. Equation (114) reduces 
to the following complex system of linear equations: 

N(N+ 1 ) 

£ Hi,j)T(j) = U(i). (118) 

3 = 1 

Since the index i is allowed to vary over the same range as j, N(N + 1) systems of equa- 
tions must be solved (in contrast to IV + 1 systems of equations for the PW method). 
Convergence of the total elastic cross sections was used to establish the number of Gaus- 
sian quadrature points used for integration. The number of Gaussian quadrature points 
for the LS3D solution method was increased to a maximum of IV = 44 points such that 
the total elastic cross sections changed less than 1% for all reactions with energies up to 
100 GeV/n. 

3.4 Localized Contributions 

It has been observed that the transition amplitude for reactions with large on-shell 
momenta — including high energy NA reactions and AA reactions at every energy — do 
not converge efficiently if the integration ranges of both momenta and polar angles are 
not restricted to regions that give significant contributions to the LS equation. The 
momenta which give non-zero contributions are estimated from the range of the optical 
potential and tend to be localized near the on-shell momentum, k. 

The range of the potential is determined from the NN transition amplitude, nuclear 
matter densities, and kinematic factors. Consider a high energy collision of a projectile 
nucleon with a nucleus, such as 4 He, which is modeled with a harmonic-well nuclear matter 
density (Gaussian). The transition amplitude and nuclear density are Gaussian functions 
of momentum transfer, |g| = k 2 + k' 2 — 2kk' cos{9). The greatest contributions occur when 
k' k at forward scattering angles. Therefore, the range of the potential is localized near 
the on-shell momentum and drops off quickly thereafter. At lower nucleon projectile 
energies, the on-shell momentum is smaller, and the range of the optical potential drops 
off more slowly. If one considers the collision of two light nuclei, then the range of 
the optical potential is further reduced because of the additional harmonic well density, 
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which is equivalent to multiplying with another q dependent Gaussian. Hence, even at 
low energies, the momenta that give non-zero optical potentials for AA reactions tend to 
be localized near the on-shell momentum. 

The range of interaction may be estimated with a high degree of certainty using a 
forward scattering angle approach. The maximum range of the potential, i? max , may be 
estimated at the forward scattering angle by moving outward from the on-shell momentum 
until the potential becomes effectively zero. As a measure of localization, a localization 
factor, fi\oci is defined, 

, max 

h-loc = r 



where k is the on-shell momentum. Small values are indicative of potentials that decay 
rapidly from the on-shell momentum. In the examples that follow, forward scattering has 
been used for the potentials, and the kinetic energy of the projectile in the laboratory 
frame has been used. 

Fig. 4 shows the real part of the optical potential for a p + 4 He reaction at 100 MeV. 
At this energy, the NR and REL potentials are nearly equal, although there is a slight 
REL shift of the maximum toward a higher on-shell momentum. The localization factor 
for this reaction is /q oc ~ 90% and is an example of a potential that is not localized near 
the on-shell momentum. 

In contrast to the low energy NA reaction, a p + 4 He reaction at 100 GeV is shown 
in Fig. 5. There is a significant shift in the REL on-shell momentum as compared to the 
NR k. It is also obvious that the maximum magnitude of the potential at k is reduced 
in the REL case. The momenta are highly localized for this reaction with a localization 
factor of fi ioc < 5%. 

Figs. 6 and 7 show the real part of the optical potential for 4 He + 20 Ne reactions at 
100 MeV/n and 500 MeV/n. At 100 MeV/n, the localization factor is /q oc « 22%, which 
is more localized than the NA reaction at the same energy. When the energy is increased 
to 500 MeV/n, the localization factor is reduced to /q oc ~ 10%. This trend shows that the 
momenta become more localized near the on-shell momenta as the energy is increased. 
In addition, localization of momenta near k occurs at lower energies per nucleon for AA 
reactions than for NA reactions. 

Figs. 8-11 illustrate the localization of momenta for equal mass systems. The real part 
of the optical potential is shown for low and high energy 4 He + 4 He reactions in Figs. 8 
and 9. Similar results are shown for 12 C + 12 C reactions in Figs. 10 and 11. As expected, 
there is a slight difference between the NR and REL potentials at low energy. Differences 
become more obvious as the maximum of the potential decreases with increasing energy. 
There is no relativistic shift of the on-shell momentum due to kinematics for equal mass 
systems. In each case, the potential dies off rapidly from the on-shell momentum, and 
localization of momentum near k is clearly observed. 

A small percentage (10 -8 %) of the maximum value of the elastic scattering optical 
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potential is used to estimate the range of the potential, i? max , and the range of momenta 
(Sr) from the on-shell value that give non-zero contributions to the LS equation, which the 
authors define as Sr = R m — k. If k — Sr > 0, then the limits of integration are changed: 
k" — >■ [k — 5n,k + 5r\. Note that the pole occurs at k" = k, but if the integration is 
handled with Sloan’s method [39] , then the principal value singularity is removed because 
the pole is placed at the center of an even number of Gaussian quadrature points. When 
k — Sr < 0, the integration proceeds in the usual way: [0, 2k] and [2k, oo). 

High energy reactions tend be forward scattering, so it is important to carefully specify 
the angular integration range such that convergence can be achieved efficiently. The 
elastic scattering optical potential was used to estimate the maximum range of the polar 
integration angle, # max , and was used for all reactions presented in this document. 

3.5 Equal Mass Kinematics 

There are several kinematic factors that may change according to the kinematics selection: 
the on-shell momentum, k] the elementary NN on-shell momentum, k e i m ; the elementary 
(NN) reduced mass, fi e i m ; the factor that transfers the transition amplitude to the scatter- 
ing amplitude, k(dk / dE)\ and the propagator. The relativistic (REL) on-shell momentum 
reduces to the non-relativistic (NR) on-shell momentum for equal mass systems; that is, 
A: Rel = fc NR = yj (mT Lab )/2, where m is the mass of the projectile or target, and T Lab is 
kinetic energy of the projectile in the lab frame. Likewise, it follows that fc^R = fcRE L . 

To examine the effect of kinematics for equal mass systems, the LS equation is written 
as an infinite series, 

T = V + VGqT = V + VGqV + VGqVGqV + .... (120) 

Using completeness of over momenta, the matrix element of the transition amplitude may 
be written, 

<k'|T|k) =(k'|Vjk) + J (k'|U|k 1 )(k 1 |G+|k 2 )(k 2 |U|k)dk 1 dk 2 (121) 

+ J (k'|U|k 1 )(k 1 |G 0 + |k 2 )(k 2 |U|k3)(k 3 |G+|k4)(k4|U|k)dk 1 dk 2 dk3dk4 + . . . . 

After operating with the propagator, (k i |Go’|k 7 -) = G$ (k, ki)S(ki — kj), the transition 
amplitude becomes, 

(k'|T|k) =(k'|U|k) + J (k / |U|k 1 )G+(fc,A; 1 )(k 1 |U|k)dk 1 (122) 

+ j (k'\V\k !)G+(fc, A: 1 )(k 1 |U|k 2 )G+(fc, fc 2 )(k 2 |U|k)dk 1 dk 2 + . . . , 


where the momenta have been relabeled for clarity. 

The only kinematic factors that appear in the potential are k e \ m and /z e i m , but, as 
previously stated, the elementary on-shell momenta are equivalent for NR and REL kine- 
matics. Next, the scattering amplitude is obtained from the transition matrix, but only 
those factors that may lead to kinematic differences between the REL and NR scattering 
amplitudes are stated explicitly: 


f(d) oc k 


dk 

dE 


£ 


1 


h'elm 


n — 1 


n 


1 

[£(*;) -£(**)] 


(123) 


where the momenta integration variables are hi and the integration symbols have been 
suppressed. The NR scattering amplitude is expressed as 


/ ne (0)oc^2 

n= 1 


n — 1 


n— 1 


4 n 

i=l 


l 

V^WY 


and the REL scattering amplitude is 


n— 1 




n = 1 


i — 1 


[1 + h{k u k)] 

r- - m 


(124) 


(125) 


where 


h(ki, k ) 


' kf + 




k 2 + m 2 ’ 


(126) 


and A is the number of nucleons in the projectile or target. It can now be seen that, for 
equal mass systems, the REL scattering amplitude reduces to the NR amplitude when 


h(ki, k) — > 1 . 


(127) 


For cases in which ki ~ k, this condition is satisfied. 

To illustrate the effect of the kinematics, the first couple of terms of the scattering 
amplitude are examined, 


/ NR (0) oc A + 


2 A 2 

k 2 - V 2 


+ .... 


(128) 


and 


/ EEL (o) cc a + + 


(129) 
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The relevant terms for comparison are 


2 


fi = 


k 2 - k" 2 


(130) 


and 


= [1 + h(k", k)) 
12 ~ k 2 - k" 2 ’ 


(131) 


where k" represents the integration variable. 

The kinematic functions, f\ and / 2 , for a 4 He + 4 He reaction at T fab = 3.3 GeV/n 
(k ~ 5 GeV) are shown in Figs. 12-14. In Fig. 12, it appears that _/j and are identical 
for all integration values, k " . However, upon closer inspection in Fig. 13, differences 
are more clearly seen when k" is farthest from the on-shell momentum. In Fig. 14, the 
black line shows the ratio / 2 // 1 , and the intersecting red, dashed lines show that the two 
functions are identical at the on-shell momentum. Fig. 12 clearly shows that the greatest 
contribution occurs at the pole. Similar results are shown in Figs. 15-17 for a 56 Fe + 56 Fe 
reaction at TL a b = 100 GeV/n. 

In this section, it has been demonstrated that the relevant kinematic factors lead to 
NR and REL scattering amplitudes which are equal at the on-shell momentum for equal 
mass systems. Although differences between the scattering amplitudes can be observed 
when integration momenta are farthest away from the on-shell momentum, the greatest 
contributions occur near the pole. As a result, very little difference between scattering 
amplitudes will be observed for equal mass systems. 


4 Results 

In the results that follow, each model uses the same set of fundamental parameterizations 
for nuclear matter densities and the NN transition amplitude. Harmonic well and two- 
parameter Fermi (Wood-Saxon) nuclear charge data are taken from references [13, 14] 
and are normalized to matter densities as described in reference [5]. When data are not 
available for the two-parameter Fermi densities, a nuclear droplet model [40] is used for 
parameter estimates. Nuclei are assumed to be near the beta stability curve. The NN 
transition amplitude used in the current work is described in Appendix A and depends 
on parameterizations of the NN cross sections, slope parameter, and real to imaginary 
ratio of the transition amplitude. The NN cross sections are taken from reference [41], 
Parmeterizations of the slope parameters and real to imaginary ratios of the data in 
reference [42] are given in Appendix B. 

In Figs. 18-24, NA and AA elastic differential cross sections are shown at energies that 
are relevant to space radiation applications, including p + 16 0, p + 56 Fe, 4 He + 4 He, 4 He 
+ 16 0, 12 C + 56 Fe, 20 Ne + 20 Ne, and 56 Fe+ 56 Fe reactions at lab projectile kinetic energies 
of 150, 500, 1000, and 20,000 MeV/n. The projectile energy of 150 MeV/n represents the 
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lower bound on model applicability. Results are indicated non-relativistic by (NR) and 
relativistic by (REL). LS3D (REL) results are given as a solid red line; a dashed, black 
line is used for eikonal results, denoted (Eik); a green square represents the PW (NR) 
results; a solid blue circle indicates LS3D (NR) results; and a violet asterisk is for PW 
(REL) results. Note that the Coulomb interaction has not been included in this analysis. 

Excellent agreement between PW and LS3D results are seen in Figs. 18-21 for each 
kinematic selection for energies greater than 150 MeV. The p + 16 0 and p + 56 Fe re- 
actions at 150 MeV/n in Figs. 18 and 19 show slight disagreements between NR PW 
and LS3D codes and eikonal results. This is likely the result of the forward scattering 
approximation used in the eikonal method, since very light projectiles may deviate from 
forward scattering at low energy. The slight disagreement between the eikonal and NR 
PW and LS3D codes is not observed for the heavier nuclei in Figs. 20 and 21, where the 
small angle scattering approximation is more appropriate. 

The next obvious feature is that of the relativistic shift observed in Figs. 18-21. The 
magnitude of the differential cross section is larger at smaller angles as compared to the 
NR cases. The effect is more pronounced at higher energies, as expected, but is also driven 
by projectile and target mass differences. A comparison of Figs. 18 and 20 shows that 
the relativistic effect is more pronounced for the p + 16 0 reaction, which has larger mass 
difference than the 4 He + 16 O system. Ultimately, the relativistic effects can be tracked 
back to kinematic differences in the relative on-shell momentum. 

Elastic scattering for equal mass systems is depicted in Figs. 22-24. Fig. 22 shows the 
elastic differential cross section for 4 He + 4 He. At the lowest energies, the eikonal model 
appears to accurately describe the scattering at small angles, but begins to deviate from 
the other models at higher angles. The most obvious difference between Fig. 22 and Figs. 
18-21 is that there is no shift in the differential cross section at higher energies in Fig. 22. 
Moreover, there is no relativistic shift observed in Figs. 23 and 24 for the 20 Ne + 20 Ne 
and 56 Fe + 56 Fe reactions. It was demonstrated in section 3.5 that the REL scattering 
amplitude is equal to the NR amplitude at the on-shell momentum for equal mass systems 
and that the largest contributions near the pole. Therefore, there is no differential cross 
section shift associated with the REL kinematics for equal mass systems. 

All of the reactions presented in Figs. 18-24 have been extended to 100 GeV/n. The 
results at this energy are consistent with the 20 GeV/n reactions and have not been 
shown here, although, as expected, there is a larger shift in the differential cross section 
for unequal mass systems. 

As an example of the LS3D method and illustration of the relativistic shift, compar- 
isons to experimental data [43-46] are performed. Fig. 25 shows the elastic differential 
cross sections of the following reactions: (a) p + 32 S at T Lab = 1 GeV [43] (b) p + 40 Ca at 
T Lab = 500 MeV [44] (c) p + 58 Ni at T Lab = 1 GeV [45] and (d) 4 He + 40 Ca at T Lab = 347 
MeV/n [46]. NR results are indicated as a solid red line, and REL results are indicated 
with a solid blue line. In each case, there is better agreement with experiment when 
relativistic kinematics are used. Since the fundamental parameterizations are based on 
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small-angle scattering data, the results are in better agreement with the measured dif- 
ferential cross section data at forward scattering angles. Also note that spin-dependence 
and medium effects have not been included, which may account for differences between 
the model and the experimental data in Fig. 25 (d). 

5 Conclusions 

The eikonal, PW, and LS3D methods have been compared for NA and AA reactions for 
reactions relevant for space radiation applications. Numerical convergence of the eikonal 
method is readily achieved when the formulas presented herein are used for light nuclei 
(A < 16) and the momentum-space representation of the optical potential is used for 
heavier nuclei. The LS formalism has an advantage over the eikonal method in that rela- 
tivistic kinematics are easily included. The PW solution method is numerically unstable 
for reactions that have large on-shell momenta, including both high energy reactions and 
relatively low energy AA reactions, due to the highly oscillatory Legendre polynomials 
needed for convergence of these systems. To circumvent this difficulty, the 3D solution 
method was implemented. Convergence of the LS3D equation can be achieved quickly 
after identifying the integration range for momenta and polar angles that give non-zero 
contributions to the LS equation. This numerical method is also useful for obtaining con- 
vergence for the partial wave analysis; however, numerical instability still exists because 
of the Legendre polynomial oscillations. 

It was shown that the NR PW and NR LS3D methods agree with the eikonal method, 
except at very low energies for projectile nucleons, where the eikonal method is not well- 
suited. As the lab energy is increased, relativistic effects are seen as a shift in differential 
cross section resonances toward higher magnitudes and lower angles. There was no ob- 
served relativistic shift of the differential cross section for equal mass systems. It was 
shown in section 3.5 that the NR and REL on-shell momentum are equal for equal mass 
systems and that the largest contributions to the LS equation occur near the pole, k. 

Based on the results presented herein, it is recommended that the LS3D method be 
used for high energy NA and AA reactions at all energies because of its rapid numerical 
convergence and stability. Extensive comparisons between models and experimental data 
are underway. The results of that research and a verification and validation analysis of 
the models will be elucidated in a future manuscript. 

6 Appendix A 

NN transition amplitude and nuclear matter densities 

In this appendix, the position space representation and Fourier transforms of the 
NN transition amplitude and nuclear matter densities are given. Usually, harmonic- well 
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nuclear matter densities are used for A < 16, and Wood-Saxon matter densities are used 
for A > 16 [5, 16]. The harmonic- well matter density in position space is [5] 


p HW (r) = (a + fir 2 ) exp 


4s 2 


where 


and 


a = 


pHW a 3 

8s 3 


1 3 y _ 37a 2 
2 8s 2 


The normalization factor is 


P = 


P r= 


Po w a 3 7 a 2 
8s 3 16s 1 ’ 

1 

7r 3 / 2 a 3 [l + 1 7 ] ’ 


(132) 

(133) 

(134) 

(135) 


where 7 and a are parameters given in references [13, 14], s 2 = a 2 / 4 — r 2 / 6, and r p is the 
proton radius [5, 16]. The Fourier transform of the harmonic-well nuclear matter density 
is given by [5] 

3 ' 21 - ~° 2s2 (136) 


p hw (?) = pr^ 3 / v[(i + ^ 7 ) - Y^ e ~ q2 

The Wood-Saxon nuclear matter density is given as [5] 


pWS( r ) = 


P^ S 


r-R 1 

1 + e 


where the normalization is 


PT = 


47 r 


R 3 + tt 2 c 2 a R 


R is the half density radius, and ca is related to the surface diffuseness, c, by [5] 

-1 


and 


2r 

Zjl p 1 

cj= vf hi 


/3 = exp 


3/3 — 1 
3-/3 


cV 3 
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(138) 


(139) 


(140) 


where r p is the proton radius. The parameters, R and c, are given in references [13, 14]. 
The Fourier transform of the Wood-Saxon matter density is [12] 


P WS (9) 


47T 

q 


■po s <t>(q), 


(141) 
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where 


0 ( g ) = ttc a R 


— cos (qR) 
sinh(gc,47r) 

O 00 

2 ca 

~ ^ 

n= 1 


7T 0,4 sin (qR) coth(gc/i7r) 
i? sinh(gc,47r) 

(— l) n nqcA 
[{qc A y + n 2 ] 2 . ’ 


(142) 


The NN transition amplitude is [5] 


^nn(’T’) = 7" exp 


25(e nn) 


(143) 


with 


T = — 



cr(eNN) 

[2vr J B(e NN )]3/2 


[«(e NN ) + i\ , 


(144) 


where cnn is the kinetic energy of the NN system in the CM frame, m p is the proton 
mass, B(e nn) is the slope parameter, ct(onn) is the NN cross section, and k(cnn) is the 
real to imaginary ratio of the NN cross section. 

The Fourier transform of the NN transition amplitude is given as [5] 


^NN(g) 


-1 1 fcg-(e N n) 
(27t) 2 ji An 


WeNN) + ile- 5 *'™* 22 ' 2 , 


(145) 


where /1 is the reduced mass of the NN system, and k is the relative momentum in the 
NN CM frame. 


7 Appendix B 

Real to imaginary ratios and slope parameters for the proton-proton (pp) and proton- 
neutron (pn) cross sections are given in this section and are parameterized to the data 
listed in reference [42], The parameterizations are a function of proton kinetic energy in 
the lab frame, e, with energy in units of GeV. 

The pp real to imaginary ratio is 


Kpp — 


Co + Ci e, 
a+ce+ee 2 +ge 3 

l+be+de 2 +fe 3 ’ 


for 0 < e < 0.2 GeV 
for e < 0.2 < 100 GeV, 


(146) 


where c 0 = 1.0736, a = -0.8370 GeV" 1 , and a = 1.3204, b = 1.0713 GeV” 1 , c = -0.9240 
GeV” 1 , d = 1.0478 GeV” 2 , e = -0.3279 GeV” 2 , / = 0.01601 GeV” 3 , and g = 0.0025 
GeV” 3 . 
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The np real to imaginary ratio is given by 


K np 


c 0 + Cie, for 0 < e < 0.15 GeV 

for e < 0.15 < 100 GeV, 


(147) 


where c 0 = 0.7763, a = 2.2886 GeV” 1 , and a = 3.2411 x 10 s , b = 1.1305 x 10 9 GeV” 1 , 
c = -7.8409 x 10 8 GeV” 1 , d = 8.0467 x 10 7 GeV” 2 . 

The slope parameter for the pp reactions is 


B 


pp 


a + ce + ee 2 
1 + be + de 2 + /e 3 


(148) 


where a = 0.4665 fm 2 , b = 0.4319 GeV 3 , c = 0.2897 fm 2 GeV 1 , d = 5.5596 xlO 3 GeV 2 , 
e = 4.5448 xlO” 3 fm 2 GeV” 2 , and f = 3.2662 xlO” 6 GeV” 3 . For the np reactions, the 


slope parameter is 


a + ce + ee 2 + ge 3 
1 + be + de 2 + /e 3 ’ 


(149) 


where a = 0.4622 fm 2 , b = -0.1776 GeV” 1 , c = -0.0139 xlO” 3 fm 2 GeV” 1 , d = 0.0207 
GeV” 2 , e = -6.5369 xlO” 3 fm 2 GeV” 2 , f = 6.9016 xlO” 3 GeV” 3 , g = 4.9074 xlO” 3 fm 2 
GeV” 3 . 


8 Appendix C 

Kinematics 

The relative on-shell momentum in the CM frame is needed as input for the LS equa- 
tion. This appendix describes the NR and REL kinematics for reactions of the form, 

A + B C + D, (150) 

where the relative momentum is written as a function of the kinetic energy of the projectile 
in the laboratory frame. 


8.1 Non-relativistic Kinematics 

With NR kinematics, the relative momentum may be expressed as 

^ _ m B pA ~ rn A PB 

r i — , 

m,A + m B 


(151) 
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where 7774 is the mass of the projectile, P4 is the momentum of the projectile, m B is the 
mass of the target, is the momentum of the target. V, = Pa cm = -Pb cm1 where pj CM 
is the momentum of particle i in the CM frame. The momentum of particle B in the lab 
frame is Ps Lab = 0, and the relative momentum can be written as 


V t 


m B 

mA + m B 


P^Lab > 


(152) 


where p.4 Lab is the momentum of particle A in the lab frame. The kinetic energy available 
in the CM frame is 


7 cm — V^ /2fx 


T, 


m B 


Lab 


itia + m B 


(153) 


where T Lab = Pn Lab /2 mA and /x = 7774777^/(7774 + m B ). If the kinetic energy of the 
projectile in the lab frame is given, then the relative on-shell momentum may be obtained 
from 


k = | Vi 



2 pm B T Lah 
m A + m B ' 


(154) 


8.2 Relativistic Kinematics 

In this section, the relativistic momentum is expressed as a function of T^b- The four- 
momentum is denoted with a superscript, //, to avoid confusion with previous discussions 
in this document. In the lab frame, the four- momenta of particles A and B are 

Pn Lab = (^ Lab ,Pn Lab ) and P Bhab — (m B: 0), (155) 

and in the CM frame, 

Pacm (Vacmi P^cm) and P Bcm (^BcmiPbcm)) (156) 

where -E) Lab is the relativistic energy for particle i in the lab frame, E icM is the relativistic 
energy for particle i in the CM frame, Pi Lab is the three-momentum for particle i in the 
lab frame, and Pi CM is the three-momentum for particle i in the CM frame. 

The Mandelstam variable, s, is an invariant dehned by the four-vector’s length, s = 
(Pa + P B ) 2 - I 11 the CM frame, s = (E A + E B , 0) 2 , and the total energy is y/s. In the lab 
frame, 


s = (^A Lab + Pb l J 2 = rn 2 A + m 2 B + 2 m B E AlMb . (157) 

The kinetic energy is found by subtracting the rest mass from the total energy, Th Lab = 

E, 4 Lab — 777-4, SO 

s = {m A + m B ) 2 + 2m B T Aljah . (158) 
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The Mandelstam variable can be used to express the total energy as a function of the 
relative three-momentum by using the CM frame, 


where "P* = Pa cm = — P-Bcm- After performing the square and regrouping terms, 
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Figure 1: Transition amplitude for AA scattering in the AA CM frame, k (k') is the 
initial(hnal) momentum of the projectile nucleus in the AA CM frame. — k (— k') is the 
initial(hnal) momentum of the target nucleus in the AA CM frame. The internal nucleon 
momenta are denoted pi and P 2 . Pp and Pp are the momenta of the core of the Ap — 1 
virtual projectile and A T — 1 virtual target nuclei, respectively. See text for explanation. 



Figure 2: Optimum factorization for AA scattering in the AA CM frame, k (k') is the 
initial(hnal) momentum of the projectile nucleus in the AA CM frame. — k (— k') is the 
initial(hnal) momentum of the target nucleus in the AA CM frame. The internal nucleon 
momenta are denoted pi and p 2 - Pp and Pp are the momenta of the core of the Ap — 1 
virtual projectile and Ap — 1 virtual target nuclei, respectively, q is the momentum 
transfer for the AA system. See text for explanation. 
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Figure 3: Illustration of the vectors used for the nucleus-nucleus optical potential. The 
distance from the center of the projectile nucleus to a nucleon in the projectile nucleus is 
rp. Likewise, vp is the distance from the center of the target nucleus to a nucleon in the 
target. The center to center distance between nuclei is r, and r NN is the distance between 
a nucleon in the projectile to a nucleon in the target. R = r + rp is the distance from the 
center of the projectile to a nucleon in the target. 
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T Lab = 100 MeV/n (p + 4 He) 



Figure 4: Real part of optical potential for a p + 4 He reaction with TLab = 100 MeV (NR 
k ~ 346 MeV and REL k « 349 MeV) at 6 cm = 0. 


T Lab = 100 GeV/n (p + 4 He) 



Figure 5: Real part of optical potential for a TLab = 100 GeV (NR k ~ 10.96 GeV and 
REL k ~ 13.63 GeV) p + 4 He reaction at 9cm = 0. 


42 


T Lab = 100 MeV/n ( 4 He + 2 °Ne) 



Figure 6: Real part of optical potential for a TLab = 100 MeV/n (NR k ~ 1.44 GeV and 
REL k « 1.46 GeV) 4 He + 20 Ne reaction at Oqm = 0. 


T Lab = 500 MeV/n ( 4 He + 2 °Ne) 



Figure 7: Real part of optical potential for a TLab = 500 MeV/n (NR k ~ 3.23 GeV and 
REL k ~ 3.39 GeV) 4 He + 20 Ne reaction at 6qm = 0. 
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T Lab = 100 MeV/n ( 4 He + 4 He) 



Figure 8: Real part of optical potential for a T La b = 100 MeV/n {k « 0.87 GeV) 4 He + 
4 He reaction at 9 cm = 0. 


T Lab = 3.3 GeV/n ( 4 He + 4 He) 



Figure 9: Real part of optical potential for a 4 He + 4 He reaction with T Lab = 3.3 GeV/n 
(k ~ 5 GeV) at 6 *cm = 0. 
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T Lab = 150 MeV/n ( 12 C + 12 C) 



Figure 10: Real part of optical potential for a 12 C + 12 C reaction at 150 MeV/n (k ~ 3.18 
GeV) with 9 cm = 0. 



k" (GeV) 


Figure 11: Real part of optical potential for a 12 C + 12 C reaction at 10 GeV/n ( k ~ 26 
GeV) with 9 C m = 0. 
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4 4 

Kinematic Functions ( He + He) 



Figure 12: Kinematic functions, /i and / 2 , for a 4 He + 4 He reaction at TLab = 3.3 GeV/n 
(k ~ 5 GeV). The integration variable is denoted k" . 


4 4 

Kinematic Functions ( He + He) 



Figure 13: Alternate view of kinematic functions, fi and / 2 , for a 4 He + 4 He reaction 
with on-shell momentum k ~ 5 GeV. The integration variable is denoted k" . 
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4 4 

Ratio of Kinematic Functions ( He + He) 



Figure 14: Ratio of kinematic functions, / 2 //i. The two functions are identical at the 
on-shell momentum, shown with intersecting red, dashed lines. The integration variable 
is denoted k" . 


Kinematic Functions ( 56 Fe + 56 Fe) 



Figure 15: Kinematic functions, /i and / 2 , for a 56 Fe + 56 Fe reaction at TLab = 100 
GeV/n ( k ~ 384 GeV). The integration variable is denoted k". 
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Kinematic Functions ( 56 Fe + 56 Fe) 



Figure 16: Alternate view of kinematic functions, f\ and fi, for a 56 Fe + 56 Fe reaction 
with on-shell momentum k ~ 384 GeV. The integration variable is denoted k" . 


Ratio of Kinematic Functions ( 56 Fe + 56 Fe) 



Figure 17: Ratio of kinematic functions, / 2 //i. The two functions are identical at the 
on-shell momentum, shown with intersecting red, dashed lines. The integration variable 
is denoted k" . 
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Figure 18: Elastic differential cross sections for p + 16 0 reactions for projectile lab ki- 
netic energies of 150, 500, 1000, 20000 MeV. Eik. represents eikonal, LS3D represents 
three-dimensional Lippmann-Schwinger, and PW represents partial wave. Non- relativistic 
results are denoted (NR) and relativistic results are denoted (REL). 
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Figure 19: Elastic differential cross sections for p + 56 Fe reactions for projectile lab 
kinetic energies of 150, 500, 1000, 20000 MeV. Eik. represents eikonal, LS3D represents 
three-dimensional Lippmann-Schwinger, and PW represents partial wave. Non- relativistic 
results are denoted (NR) and relativistic results are denoted (REL). 
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Figure 20: Elastic differential cross sections for 4 He + 16 0 reactions for projectile lab 
kinetic energies of 150, 500, 1000, 20000 MeV/n. Eik. represents eikonal, LS3D represents 
three-dimensional Lippmann-Schwinger, and PW represents partial wave. Non-relativistic 
results are denoted (NR) and relativistic results are denoted (REL). 
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0 CM ( de g rees ) 0 CM ( de S rees ) 


Figure 21: Elastic differential cross sections for 12 C + 56 Fe reactions for projectile lab 
kinetic energies of 150, 500, 1000, 20000 MeV/n. Eik. represents eikonal, LS3D represents 
three-dimensional Lippmann-Schwinger, and PW represents partial wave. Non-relativistic 
results are denoted (NR) and relativistic results are denoted (REL). 
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Figure 22: Elastic differential cross sections for 4 He + 4 He reactions for projectile lab 
kinetic energies of 150, 500, 1000, 20000 MeV/n. Eik. represents eikonal, LS3D represents 
three-dimensional Lippmann-Schwinger, and PW represents partial wave. Non-relativistic 
results are denoted (NR) and relativistic results are denoted (REL). 
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Figure 23: Elastic differential cross sections for 20 Ne + 20 Ne reactions for projectile lab 
kinetic energies of 150, 500, 1000, 20000 MeV/n. Eik. represents eikonal, LS3D represents 
three-dimensional Lippmann-Schwinger, and PW represents partial wave. Non-relativistic 
results are denoted (NR) and relativistic results are denoted (REL). 
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Figure 24: Elastic differential cross sections for 56 Fe + 56 Fe reactions for projectile lab 
kinetic energies of 150, 500, 1000, 20000 MeV/n. Eik. represents eikonal, LS3D represents 
three-dimensional Lippmann-Schwinger, and PW represents partial wave. Non-relativistic 
results are denoted (NR) and relativistic results are denoted (REL). 
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40 , 


Lab 


500 MeV (p + Ca) 




Figure 25: Elastic differential cross sections for (a) p + 32 S at T L ab = 1 GeV [43] (b) p + 
40 Ca at T Lab = 500 MeV [44] (c) p + 58 Ni at T Lab = 1 GeV [45] and (d) 4 He + 40 Ca at 
T Lab = 347 MeV/n [46], 
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